Network Behavior in Thin Film Growth Dynamics 



T. Karabacak 

Department of Applied Science, University of Arkansas at Little Rock, Little Rock, AR 72204, USA 

H. Guclu 

School of Mathematical Sciences, Rochester Institute of Technology, Rochester, NY 14623, USA 

M. Yuksel 

Department of Computer Science and Engineering, 
University of Nevada Reno, Reno, NV 89557, USA 
(Dated: May 20, 2009) 

We present a new network modeling approach for various thin film growth techniques that in- 
corporates re-emitted particles due to the non- unity sticking coefficients. We model re-emission of 
a particle from one surface site to another one as a network link, and generate a network model 
corresponding to the thin film growth. Monte Carlo simulations are used to grow films and dynam- 
ically track the trajectories of re-emitted particles. We performed simulations for normal incidence, 
oblique angle, and chemical vapor deposition (CVD) techniques. Each deposition method leads to 
a different dynamic evolution of surface morphology due to different sticking coefficients involved 
and different strength of shadowing effect originating from the obliquely incident particles. Tradi- 
tional dynamic scaling analysis on surface morphology cannot point to any universal behavior. On 
the other hand, our detailed network analysis reveals that there exist universal behaviors in degree 
distributions, weighted average degree versus degree, and distance distributions independent of the 
sticking coefficient used and sometimes even independent of the growth technique. We also ob- 
serve that network traffic during high sticking coefficient CVD and oblique angle deposition occurs 
mainly among edges of the columnar structures formed, while it is more uniform and short-range 
among hills and valleys of small sticking coefficient CVD and normal angle depositions that produce 
smoother surfaces. 



I. INTRODUCTION 

Thin film coatings have been the essential components of various devices in industries including microelectron- 
ics, optoelectronics, detectors, sensors, micro-electro-mechanical systems (MEMS), and more recently nano-electro- 
mechanical systems (NEMS). Commonly employed thin film deposition^* 2 -* 3 - techniques are thermal evaporation, sputter 
deposition, chemical vapor deposition (CVD), and oblique angle deposition. Different than others, oblique angle de- 
position techniquei£&££ & 10 1 1 1 is typically used for the growth of nanostructured arrays of rods and springs through 
a physical self-assembly process. In many applications, it is often desired to have atomically flat thin film surfaces. 
However, in almost all of the deposition techniques mentioned above, the surface morphology generates a growth 
front roughness. The formation of growth front is a complex phenomenon and very often occurs far from equilibrium. 
When atoms are deposited on a surface, atoms do not arrive at the surface at the same time uniformly across the 
surface. This random fluctuation, or noise, which is inherent in the process, may create the surface roughness. The 
noise competes with surface smoothing processes, such as surface diffusion (hopping), to form a rough morphology if 
the experiment is performed at either a sufficiently low temperature or at a high growth rate. 

Due to its intractability, a conventional statistical mechanics treatment cannot be used to describe the complex phe- 
nomenon of surface morphology formation in thin film growth. About two decades ago, a dynamic scaling approac h 12 ! 13 
was proposed to describe the morphological evolution of a growth front. Since then, numerous modeling and experi- 
mental works have been reported based on this dynamic scaling analysis-* 3 -^. On the other hand, there has been a 
significant discrepancy among the predictions of these growth models and the experimental results publishe d 15 ! 16 ! 17 . 
For example, various growth models ha ve predictions on the dynamic evolution of the root-mean-square roughness 
(RMS), 18 which is defined as w(t) = \/[h(r, t) — (h)] 2 , where h(r,t) is the height of the surface at a position r and 
time t, and (h) is the average height at the surface. In most of the growth phenomena, the RMS grows as a function 
of time in a power law form^ 2 -* 3 -^^- 9 w ~ t@ , where (3 is the "growth exponent" ranging between and 1. (3 = for 
a smooth growth front and (3 = 1 for a very rough growth front (the RMS could be as large as the film thickness). 
Figure [T] shows a collection of experimental (3 values reported in the literature^ and compares to the predictions of 
growth models. Theoretical predictions of growth models in dynamic scaling theory basically fall into two categories. 
One involves various surface smoothing effects, such as surface diffusion, which lead to (3 < Q.25," ^ 14 ! 19 The other 
category involves the shadowing effect (which originates from the preferential deposition of obliquely incident atoms on 
higher surface points and always occurs in sputtering, CVD, and oblique angle deposition) during growth which would 
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lead to (3 = l. 21 However, it can be clearly seen in Fig. [I] that experimentally reported values of growth exponent 
(3 are far from agreement with the predictions of these growth models. 20 Especially, sputtering and CVD techniques 
are observed to produce morphologies ranging from very small to very large (3 values indicating a "non-universal" 
behavior. 
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FIG. 1: (Color online) A survey of experimentally obtained values of growth exponent f3 reported in the literature^ for different 
deposition techniques is compared to the predictions of common thin film growth models in dynamic scaling theory. Root- 
mean-square roughness (RMS) grows as a function of time in a power law form, w ~ , where f3 is the "growth exponent" 
ranging between and 1. (3 = for a smooth growth front and (3 = 1 for a very rough growth front. 

Only recently, it has been recognized that in order to better explain the dynamics of surface growth one should 
take into account the effects of both "shadowing" and "re-emission" processes j 14 i 16 i 17 i 22 i 23 As illustrated in Fig. O 
during deposition, particles can approach the surface at oblique angles and be captured by higher surface points 
(hills) due to the shadowing effect. This leads to the formation of rougher surfaces with columnar structures that 
can also be engineered to form "nanostructures" under extreme shadowing conditions, as in the case of oblique angle 
deposition that can produce arrays of nanorods and nanospr ings ^aLSaS i 10 \ 1 1 In addition, depending on the detailed 
deposition process, particles can either stick to or bounce off from their impact points, which is determined by 
a sticking probability, also named "sticking coefficient" (s). Non-sticking particles are re-emitted and can arrive at 
other surface points including shadowed valleys. In other words, re-emission has a smoothening effect while shadowing 
tries to roughen the surface. Both the shadowing and re-emission effects have been proven to be dominant over the 
surface diffusion and noise, and act as the main drivers of the dynamical surface growth front The prevailing 
effects of shadowing and re-emission rely on their "non- local" character: The growth of a given surface point depends 
on the heights of near and far-away surface locations due to shadowing and existence of re-emitted particles that can 
travel over long distances. 

Figure [3] summarizes some of the experimentally measured sticking coefficient values reported in the literature 
during evaporation^ sputtering 25£&2L2&2iL2!ii21£l anc [ CVD 3 - 3 -* 3 ^ 3 -^ ££££ £8£9 g row th Q f various thin film materials. 
Names of incident atoms/molecules on the growing film are also labeled. It can be clearly seen from Fig. [3] that 
incident particles can have sticking probabilities much less than unity in many commonly used deposition systems, 
which further indicates that re-emission effects should be taken into account in attempts for a realistic thin film growth 
modeling. 

Due to the complexity of the shadowing and re-emission effects, no growth model has been developed yet within the 
framework of dynamical scaling theory that take into both these effects and still that can be analytically solved to pre- 
dict the morphological evolution of thin film or nanostructure deposition^ Only recently, shadowing and re-emission 
effects could be fully incorporated into the Monte Carlo lattice simulation approaches j 10 i 14 i 15 i 16 i 17 i 22 i 23 i 4Q i 41 i 42 i 43 These 
simulations successfully predicted the experimental results including the (3 values reported in the literature (see 
Fig. [T]). 20 However, like in experiments, (3 values from simulations ranged all the way from to 1 depending on the 
sticking coefficients used. For example, Fig. [4] shows (3 values for a Monte Carlo simulated CVD growth obtained for 
various sticking coefficient and surface diffusion values j 15 i 17 i 22 i 23 It has been observed that re-emission and shadowing 
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FIG. 2: Surface of a growing thin film under shadowing and re-emission effects is illustrated. 
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FIG. 3: (Color online) Some of the experimentally measured sticking coefficient values reported in the literature during 
evaporation, 24 sputtering,^^* 2 ^ 2 ^^^ and cyj^S^S^SLS^ growth are shown. Names of incident atoms/molecules 
on the growing film are also labeled. In same cases, depositions were done at with substrate heating at temperatures denoted 
as T s in the figure. 



effects dominated over the surface diffusion processes due to their long-range non-local character. At small sticking 
coefficients (e.g. s < 0.5) re-emission was stronger than the roughening effects of shadowing and Monte Carlo simula- 
tions produced smooth surfaces with small /3 values. At higher sticking coefficient values, shadowing effect becomes 
the dominant process and columnar rough morphologies start to form. On the other hand, like in experiments, it was 
not possible to capture a "universal" growth behavior using Monte Carlo simulation approaches, which would lead to 
dynamically common aspects of various thin film growth processes. 

Moreover, it has been very recently revealed that shadowing effect can lead to the breakdown of dynamical scaling 
theory due the formation of a mounded surface morphology. 16 ' 17 In these studies, using Monte Carlo simulations it 
has been shown that for common thin film deposition techniques, such as sputter deposition and CVD, a "mound" 
structure can be formed with a characteristic length scale that describes the separation of the mounds, or "wavelength" 
A. It has been found that the temporal evolution of A is distinctly different from that of the mound size, or the lateral 
correlation length, £. The formation of the mound structure is due to non-local growth effects, such as shadowing, that 
lead to the breakdown of the self-affinity of the morphology described by the dynamic scaling theory. The wavelength 
grows as a function of time in a power law form, A ~ t p , where p w 0.5 for a wide range of growth conditions while 
the mound size grows as £ ~ t x l z , where 1/z depends on the growth conditions. 

In brief, conventional growth models in dynamic scaling theory can not explain most of the experimental results 
reported for dynamic thin film growth; and dynamic scaling theory itself often suffers from a breakdown if shadow- 
ing effect is present, which is the case for most of the commonly used deposition techniques. On the other hand, 
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FIG. 4: (Color online) Growth exponent (3 values for a Monte Carlo simulated chemical vapor deposition (CVD) growth 
obtained for various first-impact sticking coefficient (s) and surface diffusion (D/F) values. The sticking coefficient at the 
second impact after re-emission was set to 1. Two sample surface morphologies are also included for a small s = 0.1 (left) and 
high s = 1 (right) sticking coefficient value, which leads to a smooth and rough surface topography, respectively. 



simulation techniques were not successful in revealing the possible universal behavior in various growth processes. 
Furthermore, simulations that can successfully predict the experimental results can not always be easily implemented 
by a widespread of researchers. Therefore, there is an immense need for alternative and robust modeling approaches 
for the dynamical growth of thin film surfaces that incorporates easy-to-implement analytical and/or empirical rela- 
tions which in turn can lead to universal growth behavior aspects of thin films. In this work, we explore a radically 
new "network" modeling approach for the dynamical growth of a large variety of thin film growth systems that can 
potentially capture universal properties of film growth process and at the same time not suffer from the shortcomings 
of dynamic scaling theory and Monte Carlo Simulation approaches mentioned above. 

Network modeling pervades various areas of science ranging from sociology to statistical physics or computer science, 
see Ref. 44, and the references therein. A network in terms of modeling can be defined as a set of items, referred to 
as nodes with links connecting them. Examples of real-life complex networks include the Internet, the World Wide 
Web, metabolic networks, transportation networks, social networks, etc. Recent studies show that many of these 
networks share common properties such as having a low degree of separation among the nodes modeled as Small- 
World Networks 4 ^ and having a power-law degree or connectivity distributions modeled as Scale- Free Networks. 46 



II. ORIGINS OF NETWORK BEHAVIOR DURING THIN FILM GROWTH 



Interestingly, non-local interactions among the surface points of a growing thin film that originate from shadowing 
and re-emission effects (Fig. [2]) can lead to non-random preferred trajectories of atoms/molecules before they finally 
stick and get deposited. For example, during re-emission, the path between two surface points where a particle 
bounces off from the first and head on to the second can define a "network link" between the two. If the sticking 
coefficient is small, then the particle can go through multiple re-emissions that form links among many more other 
surface points. In addition, due to the shadowing effect, higher surface points act as the locations of first-capture and 
centers for re-emitting the particles to other places. In this manner, hills on a growing film resembles to the network 
"nodes" of heavy traffic, where the traffic is composed by the amount of particles re-emitting from the nodes. In terms 
of network traffic, nodes can be classified as: source, sink, or router. So, the initial point/hill where an atom re-emits 
can correspond to a "source" in a network, and the final point where the atom sticks/settles can be thought as a 
"sink" in the network. Similarly, the intermediate reemission points/hills can be thought as the "routers". Therefore, 
a "traffic model" for the thin film growth can then be constructed by counting the number of atoms starting from a 
point on the film and ending at another point on the film. 
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III. MONTE CARLO SIMULATIONS 

Development of network models by our approach requires the track record of the trajectories of re-emitted and 
deposited atoms/particles. Since it is not possible to experimentally track the trajectories of re-emitted and de- 
posited atoms during dynamic thin film growth, we used 3D Monte Carlo simulation approaches instead, which 
were already shown to efficiently mimic the experimental processes and predict the correct dynamical growth 
morphology j 9 i 10 i 16 i 17 i 22 i 23 i 41 i 42 In these simulations, each incident particle (atom/molecule) is represented with the 
dimension of one lattice point. As substrate, we used a7Vx7V = 512x512 size lattice with continuous boundary 
conditions. A specific angular distribution for the incident flux of particles is chosen depending on the deposition 
technique being simulated. During normal angle deposition, all the particles are sent from the top along the substrate 
normal (polar angle 6 = 0°), while during oblique angle deposition simulations we used a grazing incidence flux where 
all particles are emitted at a = 85° angle from the substrate normal. For CVD, the incident flux had an angular 
spread according to the distribution function dP \0 <j)) / 'd£t = cosO/tt, where <\> is the azimuthal angle. 40 

At each simulation step, a particle is sent toward a randomly chosen lattice point on the substrate surface. Depend- 
ing on the value of sticking coefficient (s), the particle can bounce off and re-emit to other surface points. Re-emission 
direction is chosen according to a cosine distribution centered around the local surface normal/ At each impact, 
sticking coefficient can have different values represented as s n , where n is the order of re-emission (n = being for 
the first impact). In this study, we use a constant sticking coefficient value for all impacts (i.e., s n = s for all n) 
during a given simulation, which is a process also called "all-order re-emission" . 40 In all the emission and re-emission 
processes shadowing effect is included, where the particles trajectory can be cut-off by long surface features on its way 
to other surface points. After the incident particle is deposited onto the surface, it becomes a so called "adatom". 
Adatoms can hop on the surface according to some rules of energy, which is a process mimicking the surface diffusion. 
However, as noted before, non-local processes of re-emission and shadowing are generally dominant over local surface 
diffusion effects. Therefore, in this work we did not include surface diffusion in order to better distinguish the effects 
of re-emission and shadowing effects. After this deposition process, another particle is sent, and the re-emission and 
deposition are repeated in a similar way. 

In our simulations, deposition time t is represented by number of particles sent to the surface. Final simulation 
time (total number of particles sent) for all the simulations was tfi na \ = 25 x 10 7 . Because of re-emission, deposition 
rate and therefore average film thickness (d) depended on the sticking coefficient s used, and changed with simulation 
time t approximately according to d « ts/N 2 , where lattice size N was 512. 

Furthermore, in our simulations, trajectories of particles during each re-emission process can be tracked in order to 
reveal the dynamic network behavior in detail. When the simulation time reaches a pre-set value that we called the 
"snapshot state" , then we label each particle sent to the surface and start recording the coordinates of lattice point 
where the particle impacts and also the lattice point where it is re-emitted and makes another impact. Therefore, 
especially a small sticking coefficient particle can potentially make multiple re-emissions among the surface points 
and have multiple trajectory data. In order to increase the number of trajectory data for a better statistical analysis 
while keeping the surface morphology unchanged, we cancelled the final deposition of the particle sent during the 
trajectory data collection process. In other words, when the simulation time reached the pre-set value, particles were 
still being sent for re-emission and collection of trajectory data, however, they were not depositing to the surface, 
therefore not changing the surface morphology. We collected the trajectory data of about 106 re-emitted particles 
for each snapshot state. We did not include the trajectory data of particles as they re-emitted into the space or if 
they cross the lattice boundaries, since cross-boundary particles can lead to an artificially long trajectories due to 
the continuous boundary conditions used. All the simulations results are average of 10 runs (realizations), each time 
using a different seed number for the random number generator. 

IV. RESULTS AND DISCUSSIONS 

Figure [5] shows the snapshot top view images of two surfaces simulated for a CVD type of deposition, at two different 
sticking coefficients. Figure [5] also displays their corresponding particle trajectories projected on the lateral plane. 
Qualitative network behavior can easily be realized in these simulated morphologies as the trajectories of re-emitted 
atoms "link" various surface points. It can also be seen that larger sticking coefficients [Fig. [51(b) and (d)] leads to 
fewer but longer range re-emissions, which are mainly among the peaks of columnar structures. Therefore, these higher 
surface points act as the "nodes" of the system. This is due to the shadowing effect where initial particles preferentially 
head on hills. They also have less chance to arrive down to valleys because of the high sticking probabilities (see for 
example particle A illustrated in Fig.[2|). On the other hand, at lower sticking coefficients [Fig. (SJa) and (c)], particles 
now go through multiple re-emissions and can link many more surface points including the valleys that normally 
shadowed by higher surface points (particle B in Fig. [2|). This behavior is better realized in "surface-degree" and 
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their corresponding height matrix plots of Fig. [6] measured for CVD grown films at two different sticking coefficients 
s = 0.1 and s = 0.9. The high values (darker colors) in surface-degree plots correspond to the highly connected surface 
sites where these sites get or re-distribute most of the re-emitted particles. At smaller sticking coefficients [Fig.Efa)], 
which leads to a smoother morphology, surface-degree values are quite uniform indicating a uniform re-emission 
process among hill-to-hills and hill- to- valleys. On the other hand, at high sticking coefficients [Fig. E^b)], the high 
degree nodes are mainly located around the column borders suggesting a dominant column-to-column re-emission. 
This is consistent with the shadowing effect where columns capture most of the incident particles because of their 
larger heights, and also their borders are more likely to re-distribute the particles towards the neighboring column 
sides because of the re-emission process used (i.e., cosine distribution centered along the local surface normal). 




FIG. 5: Top view images from Monte Carlo simulated thin film surfaces grown under shadowing, re-emission, and noise effects 
(no surface diffusion is included in these simulations) for sticking coefficients (a) s = 0.1 and (b) s = 0.9 and with unity sticking 
coefficient at the second impacts. Each image corresponds to a 128 x 128 portion of the total lattice. The incident flux of 
particles has an angular distribution designed for chemical vapor deposition (CVD). Corresponding projected trajectories of 
the re-emitted particles are also mapped on the top view morphologies for (c) s = 0.1 and (d) s — 0.9. Qualitative network 
behavior can be seen among surface points linked by the re-emission trajectories. 

Another interesting observation revealed in our Monte Carlo simulations was the dynamic change of network 
behavior on the trajectories of re-emitted particles. Figure [7] shows top view images and their corresponding particle 
trajectories obtained from the CVD simulations for a sticking coefficient of s = 0.9, but this time at different film 
thicknesses that is proportional to the growth time. The dynamic change in the network topography can be clearly 
seen: at initial times, when the hills are smaller and more closely spaced, the re-emitted particles travel from one 
hill to another one or to a valley. However, as the film gets thicker, and some hills become higher than the others 
and get more separated, particles travel longer ranges typically among these growing hills. The shorter hills that get 
shadowed become the valleys of the system. 

It is expected that this dynamic behavior should be strongly dependent on the values of sticking coefficients and 
angular distribution of the incident flux of particles, which determine the strength of re-emission and shadowing 
effects, respectively. In other words, each deposition technique and material system can have different dynamic 
network behavior that can lead to various kinds of network systems. For example, as we will show later, the dynamic 
network among the surface points of a mounded CVD grown film can be quite different than the one among the 
nanorod and nanospring structures formed in an oblique angle deposition system, where the shadowing effect is most 
dominant, and also the one during normal angle evaporation, where shadowing effect is almost absent (re-emitted 
particles during normal angle deposition can still lead to a minimal short-range shadowing effect). 

In order to make a more quantitative analysis on the network characteristics of thin film growth dynamics, in Fig. [HI 
we plotted the degree distributions P(k) (i.e., proportional to the percentage of surface points having "degree (fc)" 
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(a)CVDs = 0.1 (b)CVDs = 0.9 



FIG. 6: (Color online) Height matrix and corresponding surface-degree values are plotted for CVD grown films with sticking 
coefficients (a) s = 0.1 and (b) s = 0.9. Total lattice size is 512 x 512 and simulation time for these snapshot states was 
t = 23.75 x 10 7 particles. 



d=7\ lattice units d = 359 lattice units d = 1 754 lattice units 




FIG. 7: First row: Top view images from Monte Carlo simulated thin film surfaces for a CVD growth with s = 0.9 at different 
film thicknesses d, which is proportional to growth time. Bottom row: Corresponding projected trajectories of the re-emitted 
particles qualitatively show the dynamic change in the network topography. 



number of links through incoming or outgoing re-emitted particles), average distance (Ik) versus degree k (i.e., the 
average "lateral" distance particles travel that are re-emitted from/to surface sites having k number of links), and 
distance distributions P(l) (i.e., proportional to the probability of a re-emitted particle traveling lateral distance I) 
for Monte Carlo simulated normal incidence evaporation, oblique angle deposition, and CVD thin film growth for 
various sticking coefficients. The left and right columns in Fig. [8] correspond to the initial (thinner films) and later 
(thick films) stages of the growth times, respectively. First, the comparison of degree distributions [Fig. [51(a) and 
(b)] of normal incidence and oblique angle growth reveals that independent of the most sticking coefficients used 
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and also their growth time, universal behavior exists for both deposition techniques: There is an exponential degree 
distribution for normal angle evaporation (confirmed in the semi- log plots, not shown here), while this behavior is 
mainly power-law for oblique angle deposition with an exponential tail. Interestingly, quantitative values of degree 
distributions for both normal and oblique angle depositions also seem to be independent of the sticking coefficient 
used, which becomes clearer at later stages of the growth [Fig. Efb)], leading to two distinct distributions for each 
deposition. The power-law observed in degree distribution of oblique angle deposition has a P(k) ~ k~ 2 behavior 
apparent at later stages. All these suggest the possibility of a universal behavior in normal and oblique angle growth 
independent of the sticking coefficient. This is quite striking since each different sticking coefficient corresponds to a 
different type of morphological growth (i.e., smoother surfaces for smaller sticking coefficients and rougher surfaces 
for higher sticking coefficients), yet the degree distribution in network traffic of re-emitted particles seems to reach a 
unique universal state. 




FIG. 8: (Color online) Behavior of degree distributions P(k) (top row), average distance (Ik) versus degree (middle row), and 
distance distributions P(l) (bottom row) for network models of a Monte Carlo simulated normal incidence evaporation (AO), 
oblique angle deposition (A85), and CVD thin film growth for various sticking coefficients s and for two different deposition 
time t (left column: t = 1.25 x 107 particles, and right column: t = 23.75 x 10 7 particles) are shown. 

As can be seen in Fig. [81(a) and (b), the re-emission process which is the dominant process in normal angle growth 
promotes an exponential degree distribution; while shadowing which is the governing effect during oblique angle 
deposition leads to a power-law distribution. On the other hand, CVD shows an exponential degree distribution at 
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initial times of the growth, while it becomes closer to power-law type for higher sticking coefficients s > 0.5. This 
is believed to be competing forces of re-emission and shadowing effects, where the re-emission is more dominant for 
smaller sticking coefficients and at initial times of the growth when the film is smoother, leading to an exponential 
degree distribution. However, shadowing effect originating from the obliquely incident particles within the angular 
distribution of CVD flux can lead to a power-law behavior at higher sticking coefficients especially when the film gets 
rougher at later stages of the growth. A power-law degree distribution corresponds to a more correlated network that 
is consistent with the long-range, column-to-column traffic observed in surface-degree plots of high sticking coefficient 
CVD above [Fig. [61(b)]. It is also realized that especially for high sticking coefficients, there exist high degree nodes 
represented with data points at the tails of the degree distributions. These relatively small percentage but highly 
connected nodes are mainly located at the column edges as seen in surface degree plot of Fig. [61(b) and are likely 
to be the "hubs" of the network. Therefore, briefly, degree distribution during CVD growth can be similar to the 
universal line of normal incidence growth for smaller sticking coefficients (s < 0.5) showing an exponential behavior 
with a short range network traffic; or it can converge to the universal power-law degree distribution of oblique angle 
deposition for higher sticking coefficients (s > 0.5) leading to a highly correlated network driven mostly at column 
edges. 

In addition, it is revealed from average distance versus degree plots of Fig. [81(c) and (d) that nodes with high degree 
are mainly linked with long-distance surface points. Independent of the deposition method, the average distance 
changes with degree k according to a power law behavior, where the value of exponent increases as the flux become 
more oblique [i.e., AO — > A85 — > CVD in Fig. (51(c) and (d)], sticking coefficient increases, and the film gets thicker 
[i.e., Fig. [51(c) and (d)]. In other words, when shadowing effect becomes more dominant and film morphology gets 
more columnar, high degree nodes can exchange atoms with longer distance surface points. This also implies that 
high degree nodes placed on column edges [Fig. [61(b)] of high sticking coefficient growth are more likely to transfer 
particles with far away other column edges as well. This process is further supported by the distance distribution plots 
of Fig. [51(e) and (f) where as the sticking coefficient is increased (i.e., less re-emission) and more obliquely deposited 
particles are introduced (more shadowing effect), more percentage particles start to travel longer distances. On the 
other hand, for smaller sticking coefficients and normal angle deposition, where the re-emission effect is more dominant, 
average distances particles travel from high degree nodes become significantly less compared to high sticking coefficient 
and oblique angle depositions. This suggests networking during re-emission dominated growth occurs mainly among 
smaller size hills and valleys, consistent with the surface-degree plot of Fig. (6](a). 

A more interesting universal behavior is observed in a "weighted and scaled average distance" versus degree plots 
of Fig. (9](a) and (b). Here we re-scale average distance for nodes with degree fc, (Ik), first with degree fc, then with the 
average distance value of all nodes (I) (average distance of all links), and plot (Ik) /(k (I)) versus k. After re-scaling, 
independent of the deposition technique used, sticking coefficients, and the growth time, all curves fall on a similar 
line obeying a power-law behavior with (Ik) / (k (I)) ~ k~ 1,2 . The origin of the -1.2 value of the exponent is not clear 
and is under investigation. 
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FIG. 9: (Color online) Weighted average distance (Ik) /(k (I)) versus degree k for network models of a Monte Carlo simulated 
normal incidence evaporation (AO), oblique angle deposition (A85), and CVD thin film growth for various sticking coefficients 
s and for two different deposition time t [(a) t = 1.25 x 10 7 particles, and (b) t = 23.75 x 10 7 particles] are shown. 



Another universal behavior is observed in distance distribution plots: Independent of sticking coefficients, normal 
incidence growth shows a power-law behavior with P(l) ~ /~ 3 . A similar power-law behavior with an exponent of -2.75 
has been observed in the distance distribution plots during a normal incidence growth simulation with re-emission 
(p. 83 of Ref.14). The authors of that work did not use a snapshot state approach, surface morphology continuously 
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changed, and therefore they measured a kind of average distance distribution of the whole growth simulation. However, 
their exponent value is still close to our results and agrees with our findings that dynamic network behavior during 
normal incidence deposition does not change significantly due to the relatively smooth morphology throughout the 
growth. On the other hand the behavior in distance distribution plots is exponential for oblique angle deposition 
(confirmed in the semi- log plots, not shown here). CVD has a power law behavior similar to that of normal incidence 
growth with P(l) ~ l~ 3 at smaller coefficients and at initial times of the growth, and becomes exponential similar to 
oblique angle deposition at higher sticking coefficients apparent especially later stages of the growth. 



V. CONCLUSION 



In conclusion, we presented a new network modeling approach for various thin film growth techniques that incor- 
porates re-emitted particles due to the non-unity sticking coefficients. We define a network link when a particle is 
re-emitted from one surface site to another one. Monte Carlo simulations are used to grow films and dynamically 
track the trajectories of re-emitted particles. We performed simulations for normal incidence, oblique angle, and CVD 
techniques. Each deposition method leads to a different dynamic evolution of surface morphology due to different 
sticking coefficients involved and different strength of shadowing effect originating from the obliquely incident par- 
ticles. Traditional dynamic scaling analysis on surface morphology cannot point to any universal behavior. On the 
other hand, our detailed network analysis reveals that there exist universal behaviors in degree distributions, weighted 
average degree versus degree, and distance distributions independent of the sticking coefficient used and sometimes 
even independent of the growth technique. We also observe that network traffic during high sticking coefficient CVD 
and oblique angle deposition occurs mainly among edges of the columnar structures formed, while it is more uniform 
and short-range among hills and valleys of small sticking coefficient CVD and normal angle depositions that produce 
smoother surfaces. 
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